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Abstract 

The screening of a Coulomb potential by superstrong magnetic 
field is studied. Its influence on the spectrum of a hydrogen atom is 
determined. 
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I will discuss recently solved Quantum Mechanical - Quantum Electrody- 
namical problem in this lectures. It was solved numerically in papers [H [2], 
and then analytical solution was found in papers |3l H]. 

We will use convenient in atomic physics Gauss units: = a = 1/137. 
Magnetic fields B > mle^ we will call strong while B > ml/e^ we will 
call superstrong. Important quantity in the problem under consideration is 
Landau radius = 1/ \/eB called magnetic length in condensed matter 
physics. 

Let us consider hydrogen atom in external homogeneous magnetic field 
B. At strong B Bohr radius is larger than an, so there are two time scales 
in the problem: fast motion in the perpendicular to magnetic field plane and 
slow motion along the magnetic field. That is why adiabatic approximation 
is applicable: averaging over fast motion we get one- dimensional motion of 
electron along the magnetic field in effective potential 
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The energy of a ground state can be estimated as 



as 




Eo = -2m\ U{z)dz\ ^ -me'^ln\B/m^e^) (2) 



and it goes to minus infinity when B goes to infinity. 

We will see that radiative corrections qualitatively change this result: 
ground state energy goes to finite value when B goes to infinity. This happens 
due to screening of Coulomb potential. 

Since at strong B reduction of the number of space dimensions occurs 
and motion takes place in one space and one time dimensions it is natural 
to begin systematic analysis from QED in D = 2. At tree level Coulomb 
potential is 

m ^ AoCk) = ^ , (3) 
while taking into account loop insertions into photon propagator we get: 
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where photon polarization operator equals: 
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fJ,U 
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t = —k"^ /Am?, and dimension of charge g equals mass in D = 2. 

Taking = (0, /cy), fc^ = — A;| for the Coulomb potential in the coordinate 
representation we get: 
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dk J 2n 



kl + Ag'^P{kyAw? 



(7) 



and the potential energy for the charges +g and —g is finally: V{z) = 
-g^{z). 

In order to perform integration in ([7]) we need simplified expression for 
P(t). Taking into account that asymptotics of P{t) are 



Pit) 



P , t < 1 
1 , t > 1 



let us take as an interpolating formula the following expression: 

2t 
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The accuracy of this approximation is better than 10%. 
Substituting ([9]) into ([7]) we get: 
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A| ^ A;J + 6m2 + 4g^ 



Akn z 



dk\\ 
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'^1^1 + n^^T^ A ^ exp(-y^6m2 + 
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In the case of heavy fermions (m ^ g) the potential is given by the tree 
level expression; the corrections are suppressed as g'^ /m?. 
In the case of light fermions {m <^ g): 

{ 7re-29l^l 2 < i In 

m<^g I -2Txg[^)\z\ , ^ » i In 

For m = we have Schwinger model - the first gauge invariant theory 
with a massive vector boson. Light fermions make a continuous transition 
from m > g to m = case. The next two figures correspond to g = 0.5, 
m = 0.1. The expression for V contains P. 
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Fig.l: Potential energy of the charges +g and —g in D = 2. The solid 
curve corresponds to P; the dashed curve corresponds to P. 
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Fig. 2: Relative difference of potential energies calculated with the exact 
and interpolating formulae for the polarization operator for g = 0.5, m = 0.1. 

To find the modification of Coulomb potential in D = 4 we need an 
expression for 11 in strong B. 

One starts from electron propagator G in strong B. Solutions of Dirac 
equation in homogenious constant in time B are known, so one can write 
spectral representation of electron Green function. Denominators contain 
k"^ — iv? — 2neB, and for B >> m^/e and A;j| << eB in sum over levels 
lowest Landau level (LLL, n = 0) dominates. In coordinate representation 
transverse part of LLL wave function is: ~ exp((— — y'^)eB) which in 
momentum representation gives ~ exp{{—k1 — ky)/eB) (gauge in which 

A = l/2[Bxr\ is used). 

Substituting electron Green functions into polarization operator we get: 



n 




X exp{ 



{q + k)l + {q + k)l 



)dqodqzjf. 
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e^B X exp{ 




)xn;.2j(A;||^fc,); 



(12) 
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m + ^1) (1 - t In (g)) + M-exp (-^j P (J, 

e^*^llMA;||d2A;^/(27r)3 



= 47re 



(13) 



/ 



+ + ^exp(-A;i/(2eS))(A;J/2m^)/(3 + A:J/2m2) ' 

(14) 
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\z\ 



(15) 



For magnetic fields B <^ 37rm^/ the potential is Coulomb up to small power 
suppressed terms: 



e^B < ml 



e 

\z\ 



1 + 



e'B 



mi 



(16) 



in full accordance with the D = 2 case, where g"^ plays the role of e^B. 
In the opposite case of superstrong magnetic fields B Sirml/e^ we get: 



r e^(-,/(i7^|.|) 1 In (.f^) > > 1 

e 



I I m 



V{z) = -e^{z) 



(18) 
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Fig. 3: Modified Coulomb potential at B = lO^^G (blue) and its long 
distance (green) and short distance (red) asympotics. 
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Fig. 4: Relative accuracy of analytical formula for modified Coulomb 

potential at B = lO^^G. 

Spectrum of Dirac equation in constant in space and time magnetic field 
is well known: 

el = ml + pI + {2n + 1 + (T,)eB , (19) 



7 



n = 0,l,2,3, cr^^il. For B > B^r = ml/e the electrons are relativistic 
with only one exception: electrons from lowest Landau level {n — 0, — 
— 1) can be nonrelativistic. 

In what follows we will study the spectrum of electrons from LLL in the 
Coulomb field of the proton modified by the superstrong B. 

Spectrum of Schrodinger equation in cylindrical coordinates (p, z) in the 
gauge, where A — |[-Bf] is: 



\m\+'m + l + az\eB ^ pi 
2 I nif, 2mp 



(20) 



LLL corresponds to Up = 0, az = —1, m = 0, —1, —2, .... 

A wave function factorizes on those describing free motion along a mag- 
netic field with momentum and those describing motion in perpendicular 
to magnetic field plane: 



Rom{p)= vr(2al,)i+l™l(|m|!) 



-1/2 



\m\ {imip-p^/{4a%)) 



(21) 



Now we should take into account electric potential of atomic nuclei situ- 
ated aX p = z = 0. For an clb adiabatic approximation is applicable and 
the wave function in the following form should be looked for: 



*nOm-l = Rom{p)Xn{z) , 



(22) 



where Xn{z) is the solution of the Schrodinger equation for electron motion 
along a magnetic field: 



1 



+ Ueff{z) 



Xniz) = EnXn{z) 



(23) 



2m dz"^ 

Without screening the effective potential is given by the following formula: 



eff\ 



d'^p , 



For 1 2; I » an the effective potential equals Coulomb: 

z an 

and effective potential is regular at 2; = 0: 



(24) 
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(25) 
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an 



(26) 



Since Ueff{z) = Ueff{—z), the wave functions are odd or even under 
reflection z — )■ —z; the ground states (for m = 0, —1, —2, ...) are described by 
even wave functions. 

To calculate the ground state of hydrogen atom in the textbook "Quan- 
tum Mechanics" by L.D. Landau and E.M.Lifshitz the shallow-well approxi- 
mation is used: 



E'^ = -2m, 



as 



U{z)dz 



= -{meey2)ln''{B/{my)) 



(27) 



Let us derive this formula. The starting point is one- dimensional Schrodinger 
equation: 

1 d"^ 

;X{z) + U{z)x{z) = Eoxiz) 



2/i dz2' 

Neglecting Eg in comparison with U and integrating we get: 



a 

X'{a) = S/^y^ U{x)x{x)dx , 



(28) 



(29) 



where we assume U{x) = U{—x), that is why x is even. 

The next assumptions are: 1. the finite range of the potential energy: 
U{x) 7^ for a > a; > —a; 2. x undergoes very small variations inside the 

well. Since outside the well x(a^) ~ e 

-y/MEo\ readily obtain: 



\En 



2ii 



U {x)dx 



For 



Upla^ -C 1 



(30) 



(31) 



(condition for the potential to form a shallow well) we get that, indeed, 
\Eo\ <^ \U\ and that the variation of x inside the well is small, Ax/x ~ 
li\U\a^ <S 1. Concerning the one-dimensional Coulomb potential, it satisfies 
this condition only for a <C l/(mee^) = ob- 

This explains why the accuracy of log^ formula is very poor. 



9 



Much more accurate equation for atomic energies in strong magnetic field 
was derived by B.M.Karnakov and V.S.Popov [S]. It provides a several per- 
cent accuracy for the energies of EVEN states for H > 10^ {H = B/{mle^)). 

Main idea is to integrate Shrodinger equation with effective potential from 
a; = till X = z, where an « z « clb and to equate obtained expression 
for 'x'{z)/x{z) to the logarithmic derivative of Whittaker function - the so- 
lution of Shrodinger equation with Coulomb potential, which exponentially 
decreases at z » as- 



2 In — + \n2 - + \m\) + 0{aH/z) = 
\(if[ / 

21n(^^^ + X + 2\nX + 2ilj(^l-^ + 47 + 2 In 2 + 0(^/aij) ,(32) 



E = -(meeV2)A' 
The energies of the ODD states are: 



odd 



' 2n2 



+ 



n = 1,2, 



(33) 



(34) 



So, for superstrong magnetic fields B ~ m^/e^ the deviations of odd states 
energies from the Balmer series are negligible. 

When screening is taken into account an expression for effective potential 
transforms into 



TT r \ 2 / |-R0m(p)P ,2 



l-e 



(35) 



y/p^Tz^ 

The original KP equation for LLL splitting by Coulomb potential is: 

In(iJ) = A + 21nA + 2^/' (^1 - +ln2 + 47 + V(l + |m|) , (36) 

where ^/'(x) is the logarithmic derivative of the gamma function; it has simple 
poles at x = 0, —1, —2, .... 

The modified KP equation, which takes screening into account looks like: 



In 



H 



\ 371 / 



A + 21nA + 2^ (^1 - +ln2 + 47 + ^/'(l + |m|) , (37) 
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E = — (mee^/2)A^. In particular, for a ground state A = —11.2, Eq = —1.7 
keV. 
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Fig.5: Spectrum of hydrogen levels in the limit of infinite magnetic field. 
Energies are given in rydberg units, Ry = 13.6 eV. 



In conclusion, 

1. analytical expression for charged particle electric potential in c? = 1 is 
given; for m < g screening take place at all distances; 

2. analytical expression for charged particle electric potential ^{z,p = 0) 
at superstrong B a.t d = 3 is found; screening take place at distances \z\ < 

l/lTLe] 

3. an algebraic formula for the energy levels of a hydrogen atom originat- 
ing from the lowest Landau level in superstrong B has been obtained. 
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tality. I was supported by the grant RFBR 11-02-00441. 
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